Global generalized synchronization in networks of different time-delay systems 
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We show that global generalized synchronization (GS) exists in structurally different time-delay 
systems with different fractal dimensions, which emerges via partial GS in symmetrically coupled 
network with ring, array, global and star configurations. We find that there exists a smooth trans- 
formation even for systems of different types, which maps them to a common GS manifold. The 
asymptotic stability of all the systems in the network guarantees the existence of GS and we confirm 
it by using their maximal transverse Lyapunov exponent (MTLE). In addition, an analytical stabil- 
ity condition using the Krasvoskii-Lyapunov theory is deduced. Understanding the onset of GS in 
such systems may lead to an effective control of pathological synchronous activity with extremely 
small perturbations to appropriate variables. 

PACS numbers: 05.45.Xt,05.45.Pq,02.30.ks,89.75.-k 



Synchronization is a ubiquitous nonlinear phenomenon 
serving as a platform for information processing in di- 
verse natural and man made systems [l|, 0| ■ It has been 
investigated mainly in identical systems and systems with 
parameter mismatch, with rare exceptions of distinctly 
(structurally) nonidentical systems. However, in reality 
very often synchronization emerges in distinctly noniden- 
tical systems such as respiratory arrhythmia between car- 
diovascular and respiratory systems [3| , visual and motor 
systems [1] , paced maternal breathing on fetal , differ- 
ent populations of species 0,0], m epidemics jgj, in cli- 
matology and many more. Considering the coherent 
coordination of living systems involving multiple organs 
such as brain, heart, lungs, limbs, etc., or machines con- 
sisting of distinct parts, cooperative evolution of distinct 
and often time-delayed systems is essential and challeng- 
ing. 

Among various kinds of synchronization admitted by 
coupled nonlinear dynamical systems [l|, 0, the intri- 
cate phenomenon of generalized synchronization (GS) 
refers to (static) functional relationship between inter- 
acting systems |UHLa |. While the phenomenon of GS 
has been well understood in unidirectionally coupled sys- 
tems [Iol - [l3j . it remain still in its infancy in bidirection- 
ally coupled systems and in particular there exists only 
very limited results on GS, even in systems with param- 
eter mismatches 

EMI, 

and particularly in structurally 
different (nonidentical) systems with different fractal di- 
Thus, in general, the notion of GS in 
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mensions 

mutually coupled systems needs much indepth investiga- 
tion and in particular in distinctly different systems with 
different fractal dimensions including delay systems. In- 
deed, recent investigations have revealed that GS is es- 
sentially more likely to occur in complex networks (even 
with identical nodes) [l5| due to the large heterogeneity 
(degree distribution) of many natural networks }16] . 



It is important to recall that the above mentioned 
studies [J-Q have demonstrated only PS among such 
distinctly different complex systems, while the natural 
choice of GS in such systems has been largely neglected. 
Further, depending on the relation between PS and GS in 
systems with different fractal dimensions, which remains 
unclear, our understanding on their evolutionary mech- 
anism, dynamical and functional behavior may need to 
be reinvestigated. Controlling pathological synchronous 
activity and inducing coherent coordination in paralysed 
systems may be effectively done upon understanding the 
emergence of a common GS manifold in such systems. 

The main reason to consider time-delay systems in 
this work is due to the ubiquity of time-delay systems 
as several real systems in ecology, epidemics, physiology, 
economics, engineering and control systems, etc. Q in- 
evitably require intrinsic delay (different from connection 
delay between different units) for a complete description. 
The flexibility of choosing time-delay systems with dif- 
ferent fractal dimensions just by adjusting their intrinsic 
delay times is quite attractive from a modelling point of 
view. In such a situation the most crucial aspect is to ex- 
amine the possibility (whether functional relation exists) 
to map distinctly nonidentical time-delay systems with 
different fractal dimensions to a common synchroniza- 
tion manifold in spite of their high-dimensional complex 
phase space with multiple positive Lyapunov exponents. 

In line with the above discussions, in this Letter we will 
provide a substantial extension of the formulation of GS 
to mutually coupled distinctly different delay systems, 
which also holds for systems without delay. Based on this 
we will specifically demonstrate the existence of global 
GS via partial GS in symmetrically coupled networks, 
which even consist of distinctly different time-delay sys- 
tems (Mackey-Glass [2(|, piecewise linear [23j|, threshold 
nonlinearity [13] and Ikeda ; 21j]). It is surprising that 
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FIG. 1: (color on-line) Schematic diagram of the auxiliary 
system approach for mutually coupled systems in a ring con- 
figuration for (a) N = 2 and (b) N = 4. 

there exists a common GS manifold even in an ensemble 
of distinctly different time-delay systems to achieve global 
GS in four different network (array, ring, star and all-to- 
all) configurations, i.e, there exists a function (smooth 
map) for each system, even with different fractal dimen- 
sions, in a network which maps them to the common 
GS manifold. We calculate the maximal transverse Lya- 
punov exponent (MTLE) to evaluate the asymptotic sta- 
bility of the complete synchronization (CS) manifold of 
each of the main systems with their corresponding aux- 
iliary systems starting from different initial conditions, 
which in turn asserts the stability of the GS manifold be- 
tween the main systems. Further, we will also estimate 
cross correlation (CC) and correlation of probability of 
recurrence (CPR) [22J to establish relations between GS 
and PS. 

Kocarev and Parlitz formulated that GS in drive x 
and response y configuration occurs only if the response 
is asymptotically stable, i.e. V yi(0),x(0) in the basin 
of synchronization manifold lim^oo ||y(t, x(0), yi(0)) — 
y(t,x(0),y2(0))|| = 0, a mathematical formulation of the 
concept of auxiliary system approach jllfl . Now we give a 
substantial extension of this formulation to mutually cou- 
pled different dynamical systems with delay represented 
as 

x = f(x,x r ,u), y = g(y,yr,v), f^g (l) 

where x, x T G R", y, y r £ R m , r € R and u, v e R fc , k < 
to, n. Ui = —Vi = hi (x(t, xo), y(£, yo)) correspond to the 
driving signals. System ([1]) is in GS if there exists a trans- 
formation H : R n , R m -> C R™ +m . This implies that 
the synchronization manifold M = {(x, y) : H(x, y) = 0} 
is such that V x(0),y(0) which lies within the subset 
B = B x x B y C R" x R m approaches M C B, so that 
M is an attracting invariant manifold. 

Thus, GS exists in system ([1]) only if 
both coupled systems are asymptotically sta- 
ble such that V(y<(0),Xf(0)) C B, i = 1,2, 
lim t _> oo ||y(f,xi(0),y 1 (0)) - y(t,x 1 (0),y 2 (0))|| = 
and lim^oo ||x(t,xi(0),yi(0)) - x(i, x 2 (0), yi (0))|| = 0. 
This is a mathematical formulation of the auxiliary 
system approach to system (JIJ, whose schematic dia- 
grams with the main and auxiliary systems in a ring 
configuration are depicted in Figs. [TJ, and [Tb, respec- 
tively. Therefore trajectories of ([!]) starting from B 
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FIG. 2: (color on-line) (a) MTLEs and (b) CC, CPR of the 
main and auxiliary systems of coupled MG-PWL systems. 

asymptotically reach M defined by the transformation 
function H(x, y), which can be smooth if the systems 
((T|) uniformly converge (otherwise nonsmooth) , i.e. their 
local Lyapunov exponents (LLEs) are always negative, 
to M 0. 

Now, we will demonstrate the existence of GS in sym- 
metrically coupled arbitrary network of distinctly non- 
identical time-delay systems with different fractal dimen- 
sions using the above formalism. The dynamics of the zth 
node in the network is represented as 

JV 

±i = -OiXi(t) + /3ifi(xi(t - n)) - e}^Gjj-Kj, (2) 

where i = 1,...,N, N is the number of nodes in the 
network, and /3i are constant parameters, Tj are the 
delay times, fi is the nonlinear vector function of ith 
node and G is an Laplacian, ensuring the invariant man- 
ifold M, determining the topology of the arbitrary net- 
work. To determine the asymptotic stability of each of 
nodes in the network, we define an idcntial (auxiliary) 
network (starting from different initial conditions) with 
node i represented as x'i — — ajX-(i) + /3j£-(x^(i — Tj)) — 

e^jli < -^ij(%( x i — x j) + x j) ( see Fig- [J). Similar anal- 
ysis has been carried out for N — 3, 7 and 10, to confirm 
the generality of the phenomena described in the follow- 
ing. For brevity, we are not including those cases here. 

First, we consider N — 2 mutually coupled dis- 
tinctly nonidentical time-delay systems with the nonlin- 
ear function f\{x) = (i^x^t-T^) 10 )) °^ tne Mackey-Glass 
(MG) [13] system and the piecewise linear (PWL) func- 
tion f 2 {x) = Oifx < 4/3,-1.5a:-2if-4/3 < x < -0.8, x 
if -0.8 < x < 0.8, -1.5a: + 2 if -0.8 < x < 4/3, and 
if x > 4/3 of the PWL time-delay system (23[, respec- 
tively. The parameters of the MG systems are chosen 
as ot\ = 0.5, j3\ = 1.0, and t\ — 8.5 and for the PWL 
systems we choose ai = 1.0, — L2 and r 2 = 10.0. 
For these parameter values both systems are known to 
exhibit hyperchaotic attractors with two {Dry — 2.957) 
and three (B>ky = 4.414) positive Lyapunov exponents 
(LEs) [2j, respectively. 
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At a first sight, one may think that mutually interact- 
ing systems compromise to reach a common synchroniza- 
tion manifold simultaneously for a certain critical value of 
e. But because of the distinctly different coupled systems 
with different fractal dimensions, one of the systems first 
reaches a GS manifold for a lower value of e, while the 
other system remains in a desynchronized state, which 
we refer to as a partial GS. With further increase of e the 
other system also converges at a different critical e to 
the common GS manifold achieving a global GS. Further, 
one may also expect that systems with lower dynamical 
complexity will converge to GS manifold first and then 
the higher one in the order of their degree of complex- 
ity (here in terms of the number of positive LEs). On 
the contrary, the PWL system with three positive LEs 
reaches the GS manifold first at e « 0.26 as indicated 

(2) 

by the changes in the sign of the ^mtle between 2 and 
2' in Fig. [5^ and then the MG system with only two 
positive LEs smoothly converges to the GS manifold at 



further larger e « 0.5, where A 



(i) 

MTLE 



< as shown in 



Fig. [5^ confirming the emergence of a global GS in dis- 
tinctly nonidentical time-delay systems. 

The emergence and the transition from par- 
tial to global GS is also characterized by the 
correlation coefficient (CC) given by d f = 

((x t (t)-(x z (t))){x[(t)-{x[(t)))) 



where the 



brackets 



V((^(o-<xi(t)» 2 >(K(t)-K(t)» 2 >' 

indicate time average. If the two systems are in CS state 
the CC — I, otherwise CC < 1. Further, the existence of 
PS in highly non-phase-coherent hyperchaotic attractors 
of time-delay systems are characterized by the value 
of the correlation of probability of recurrence (CPR 

« 1) mm. 

In the absence of coupling (e = 0.0), CC\^y and CCi^i 
and CPRi^v and CPR 2 ^ d2b) are nearly zero and both 

^mtle an d ^mtle > 0, indicating the absence of CS 
(GS) between the main and auxiliary systems. If we in- 
crease the coupling strength, CC^' and CPi?2,2' start 
to increase towards unity and at Ec ~ 0.26, CC'2.2' = 1 

(2) 

(CPi?2,2' = 1), where A\/ TLB < , which confirms the 
onset of GS (PS) in the PWL system while the MG 
system continues to remain in a desynchronized state 



(CC 2 



0.2 and A 



(i) 

MTLE 



> 0). Further, if we increase 



(2) 



0.5, a 



the coupling strength to a threshold value e 
global GS occurs where both CCi } y and CC-z,^ become 

unity and while both ^mtle an( ^ ^mtle ^ 0. It is also 
to be noted that both CPRi^y and CPR2.2' are also in 
agreement with their CC confirming the existence of GS 
and PS together. 

We have also analytically investigated the existence 
of partial and global GS using the Krasovskii-Lyapunov 
theory. For this purpose, we consider the difference in 
the state variables of the main and auxiliary systems 
(Aj = Xi — x'j) for small values of A*. The evolution 
equation for the CS manifold for the Mackey-Glass time- 
delay systems (Ai = x\ — x[) can be written as Ai = 
-(ai + e)Ai + /3i[/{(a; T1 )]Ai Tl , where x Tl = x(t - n) 



FIG. 3: (color on-line) (a) and (c) MTLEs and (b) and (d) 
CC, CPR of the main and auxiliary systems of four distinctly 
nonidentical time-delay systems with ring and array configu- 
rations, respectively, as a function of e. 



and Ai Tl = Ai(t — T\). The CS manifold is locally at- 
tracting if the origin of Ai is stable. Then, a continuous 
positive definite Lyapunov functional can be defined in 
the form V(t) = ~Af + fi J°_ r A? (t + 0)d6, 7(0) = 0, 
where /1 is an arbitrary positive parameter (/1 > 0). The 
Lyapunov function V(t) approaches zero as Ai — > 0. The 
derivative of V(t) along the trajectory of the CS manifold 
should be negative for the stability of the CS manifold 
Ai = 0. This requirement results in a sufficiency condi- 
tion for the asymptotic stability as (ai+e) > \/3if[(x Tl )\. 
From the form of the nonlinear function f(x) for the 
MG system, the stability condition becomes (cvi + e) > 



(l+xi 



-cx T 

J*- 



From the hyperchaotic attractor of the Mackey-Glass 

system for the above chosen parameter values one can 

find numerically the maximum value of x max ss 1.24. So 

the stability condition becomes e > \flif[(x T )\ — a sa 0.41 

(which is a sufficiency condition for global GS). From 

Fig. [51 the threshold value of the coupling strength to 

(2) 

attain GS in the Mackey-Glass systems is £c ~ 0.45 (> 
0.41) for which the stability condition is satisfied. With 
a similar procedure, the stability condition for the PWL 
systems become (a% + e) > \02f2i x n)\- N° w from the 
form of the piecewise linear function /b^c), we have 
fi( x v>) = 1-5 for 0.8 < \x\ < § and 1.0 for |x| < 0.8. Con- 
sequently, the stability condition becomes (a 2 +e) > \f3 2 \, 
corresponding to the inner regime of the nonlinear func- 
tion where most of the system dynamics is confined, for 
the asymptotic CS state A2 = 0. From Fig. [21 we find 
that CS between the piecewise linear systems {X2 and 
x' 2 ) occur for the coupling strength Ec > 0.27, which 
satisfies the sufficient stability condition for partial GS 
e > |/3 2 | -0L2 = 0.2 (0.27 > 0.2). 

Next, we demonstrate the existence of GS in four 
(only in N — 4 for clear visibility of figures depict- 
ing synchronization transitions) mutually coupled dis- 
tinctly nonidentical time-delay systems in a ring, an ar- 
ray, a global and a star configuration. In addition to the 
above two time-delay systems, now we have considered 
the third system as a threshold piecewise linear time- 
delay (TPWL) system [24] with the nonlinear function 
/ 3 (x) = AF* - Bx, where F* = — x* for x < -x*, x for 
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FIG. 4: (color on-line) (a) and (c) MTLEs and (b) and (d) 
CC, CPR of the main and auxiliary systems of four distinctly 
nonidentical time-delay systems with global and star config- 
urations, respectively, as a function of e. 

—x* < x < x* and x* for x > x*. The system param- 
eters are fixed as follows: = 1.0, f3^ = 1.2, T3 = 7.0, 
A = 5.2, B = 3.5 and x* = 0.7 and the system ex- 
hibits hyperchaotic attractor with four positive LEs. As 
a fourth system, we take the Ikeda time-delay system [2l[ 
with the nonlinear function fi(x) — sin(x(t — T4)). This 
system exhibits hyperchaotic attractor with five posi- 
tive LEs [2[ for the parameters a a — 1.0, /?4 = 5.0 and 
T4 = 7.0, which is used in the remaining part of the paper. 
CC, CPR and MTLE for mutually coupled ring of the 
above four distinctly nonidentical systems are shown in 
Figs. [5^, and [5b. ^} TLE > and low values of CCi^ and 
CPRi.i' for e = indicates the main and their auxiliary 
systems evolve independently. Increasing e results in de- 
creasing )^mtle an d increasing CCi^> and CPRi^i. The 
PWL systems with 3 positive LEs become synchronized 

(2) (2) 

first at £c = 0.088 as evidenced by ^ MTLE < 0, indi- 
cating the onset partial GS, and then the MTLE of MG 
with 2 positive LEs becomes negative {^m TLE < 0) at 

ei 1 ^ ~ 0.092 while the other two systems are not yet syn- 
chronized. Further, increase in e leads to CS of Ikeda sys- 
tems with 5 positive LEs at ei^ ~ 0.65 with \^j TLE < 
and finally the TPWL system becomes synchronized at 
ef^ w 0.75, confirming the existence of a global GS be- 
tween the four mutually coupled distinctly nonidentical 
time-delay systems in a ring configuration. In addition, 
all CCi t i> and CPRi : i> reach unity (Fig. [3b) exactly at 

the threshold values e c of ^mtle corroborating the si- 
multaneous existence of GS and PS. 

Next, we illustrate the transition from partial GS to 
global GS in a linear array of N = 4 (MG, PWL, 
TWPL and Ikeda) systems. MTLE, CC and CPR of 
all the four main and their auxiliary systems are shown 
in Figs. [3J; and [3J1. Again, all the four systems reach 
their CS (GS) manifolds for the threshold values Sc ~ 
0.09, 0.04, 0.54, 1.2, respectively, as indicated by the tran- 
sition of their \j$tle below zero. Further, GS and PS 
also occur together as indicated by CC^y and CPRi^ 
(Fig. [3ji) of the systems in the array. 

Maximal TLE's and CC and CPR of all the four sys- 
tems with global and star configurations are depicted in 



Figs. Hk, Hp and St, HH, respectively. As in the case 
of ring configurations, ^jJ TLE transits first from positive 
to negative values elucidating CS between PWL and its 
auxiliary system thereby indicating the onset of partial 
GS. This is followed by the MG systems, then by Ikeda 
systems and finally the TPWL systems for appropriate e 
confirming the existence of GS in four distinctly noniden- 
tical time-delay systems with different fractal dimensions 
in both global and star configurations. 

Further, we have also confirmed that there exists tran- 
sition from a partial to global GS in other permutations 
on the order of systems between MG, PWL, TPWL, 
Ikeda systems in array and star configurations. Follow- 
ing a similar stability analysis as above for the four cou- 
pled time-delay systems with their auxiliary systems for 
CS, one can also arrive at a sufficiency stability condi- 
tions for GS in all the above four configurations. All the 
above results have also been confirmed in a larger en- 
semble (N = 7, 10) of time-delay systems with distinct 
fractal dimensions in all the four configurations. 

In summary, we have found that there exist smooth (as 
evidenced by the smooth convergence of MTLE) trans- 
formation functions even for distinctly nonidentical time- 
delay systems with different fractal dimensions in sym- 
metrically coupled network (ring, array, all-to-all and 
star), which map each of them to a common synchro- 
nization (GS) manifold. We have also shown that the 
asymptotic stability of each of the system in the network 
guarantees the existence of GS, which is confirmed by 
their MTLEs and analytical stability conditions. In ad- 
dition, we find that GS always coexists with PS or vice 
versa in these systems using CC and CPR and hence our 
understanding on the evolutionary mechanism including 
their dynamical and functional behavior of systems with 
different fractal dimensions [3-9] should be revisited and 
improved. For instance, GS leading to pathological dis- 
orders may be effectively controlled with extremely small 
external perturbations to appropriate variables (to elec- 
trodes) without harming the subject. In this connection, 
investigations on desynchronization mechanisms and con- 
ditions in ensembles of such systems will also be crucial. 
The concept of equivalence and predictability, which is 
well understood in unidirectionally coupled nonidentical 
systems [Io| , remains an open problem for bidirectionally 
coupled systems. Further, one can extend the present 
work to complex networks and also include heterogenity 
by weighted couplings, etc., to generalize much more. 
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